Dynamic optical lattices of sub-wavelength spacing for ultracold atoms 
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We propose a scheme to realize lattice potentials of sub-wavelength spacing for ultracold atoms. 

It is based on spin-dependent optical lattices with a time-periodic modulation. We show that the 
atomic motion is well described by the combined action of an effective, time-independent, lattice 
of small spacing, together with a micro-motion associated with the time-modulation. A numerical 
simulation shows that an atomic gas can be adiabatically loaded into the effective lattice ground 
state, for timescales comparable to the ones required for adiabatic loading of standard optical lattices. 
We generalize our scheme to a two-dimensional geometry, leading to Bloch bands with non-zero 
Chern numbers. The realization of lattices of sub-wavelength spacing allows for the enhancement 
of energy scales, which could facilitate the achievment of strongly-correlated (topological) states. 


Optical lattices have allowed experiments on ultracold 
atomic gases to investigate a large range of lattice models 
of quantum many-body physics [T]. Their development 
led to the realization of strongly-correlated states of mat¬ 
ter, such as bosonic and fermionic Mott insulators, and 
low-dimensional gases [2] . In its simplest form, an optical 
lattice consists of the optical dipole potential associated 
with a standing wave of retro-reflected laser light. It can 
be described as a periodic potential V(x) = Uq cos 2 (kx), 
of spatial period d = A/2, where A is the laser wavelength 
and k = 2 tt/X. More complex optical lattices, such as su¬ 
perlattices pan or two-dimensional honeycomb lattices 
[3 El, can be generated with suitable laser configurations. 
The recoil energy E r = h 2 /(8md 2 ), where h is Planck’s 
constant and m is the atom mass, sets the natural energy 
scale for elementary processes, such as atom tunneling 
between neighboring lattice sites, as well as the temper¬ 
ature range T < E T /kB ~ 100 nK, typically required for 
quantum degeneracy. 

For a large class of models, the physical behavior is 
dictated by processes associated with even much smaller 
energies, such as super-exchange or magnetic dipole in¬ 
teractions TJ. The associated temperature scales remain 
out of reach in current experiments. In order to circum¬ 
vent this limitation, it is desirable to find novel schemes 
for generating optical lattices with spacing d e g -C A, 
in order to enhance the associated energy scale E^ s = 
h 2 /(8md 2 eS ) 0. Schemes have been proposed to gener¬ 
ate lattices of sub-wavelength spacing, based on multi¬ 
photon optical transitions |5] or on adiabatic dressing of 
state-dependent optical lattices [7] ; the realization of lat¬ 
tices with spacing d e g = A/4 was reported in Ref. [9]. An 
interesting alternative would be to trap atomic gases in 
the electromagnetic fields of nano-structured condensed- 
matter systems |TUHT2] . 

In this letter, we propose a novel scheme leading to 
lattices of spacing d e g = d/N, N being an arbitrary inte¬ 
ger, based on spin-dependent lattices with time-periodic 
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FIG. 1. Stroboscopic scheme for engineering short-spacing 
lattices, illustrated on the case N = 4. We make use of a 
periodic potential V ( x , t) of spatial period d, that is shifted of 
the distance d/N after every time step of duration T/N (blue 
curves). The effective potential V e g{x) (red curve), resulting 
from time averaging, exhibits a spatial period d e g = d/N. 


modulation. In the regime of large modulation frequency 
nMm, the atom dynamics is governed by an effec¬ 
tive static periodic potential of spacing d e g, with an 
additional micro-motion. This description is confirmed 
by a numerical simulation, which shows the possibility 
to load adiabatically the ground state of the effective 
lattice and to perform Bloch oscillations. We discuss 
the extension of the scheme to two-dimensional lattices 
with non-trivial topology. Lattices with artificial mag¬ 
netic fields, generally leading to topological bands, were 
recently realized in experiments, with standard lattice 
spacing m- For those systems, increasing the energy 
scale using short-spacing lattices could prove important 
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for creating strongly-correlated states such as fractional 
Chern insulators [18j [19]. 

A basic scheme of our method is pictured in Fig.|T| 
Consider a periodic potential V(x) of period d, which 
is abruptly shifted by the distance d/N at stroboscopic 
times t n = ( n/N)T , n £ Z, leading to a time-periodic 
potential V ( x , t) of period T. Provided that T is much 
smaller than typical timescales of atomic motion, the 
atoms experience an effective time-averaged potential 
V e g(x) = fV(x, t)dt/T. A simple calculation shows that 
V e g{x) is given by the sum of all harmonics of the po¬ 
tential V ( x ), whose orders are multiples of N [2D] . The 
effective potential V e g(x) is thus spatially periodic, of 
spatial period d e g = d/N. 

Conventional optical lattices present a spatial modu¬ 
lation proportional to the intensity pattern of interfering 
light waves, which exhibit spatial frequencies of at most 
twice the light momentum k. Thus, applying the stro¬ 
boscopic scheme in Fig. [l] to these potentials could not 
lead to effective lattices of period d e g < A/2. This re¬ 
striction does not apply to spinful particles subjected 
to spin-dependent optical lattices. As an illustration, 
consider a spin-1/2 particle evolving in the potential 
V(x) = VlCOs(2 kx)a z + Vb <t x , where a u (u = x,y,z) 
are the Pauli matrices. In a dressed state picture, the 
atom may follow adiabatically the state of lowest energy 
V.(X) = -y/V* cos 2 (2 kx) + Vq . As this potential ex¬ 
hibits harmonics of the spatial frequency 2 k of all orders, 
the lattice spacings achievable by applying the strobo¬ 
scopic scheme to V-(x) can be made arbitrarily small. 

We describe in the following a modified, more practical, 
version of this scheme, which consists of a spin-dependent 
optical lattice with smooth temporal variations, given by 
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FIG. 2. Band structure of a dynamic optical lattice of spac¬ 
ing d e g = d/4, corresponding to the parameters N = 4, and 
Vl = Vb = hQ = 200.E®®. In (a), we make use of the spa¬ 
tial and temporal translational symmetries T x , Tt and label 
the eigenstates by their quasi-momentum —k<q<k and 
quasi-energy —hfi/2 < hio < hfi/ 2. The Bloch-Floquet bands 
can be unfolded using the additional symmetry T*, leading 
to the band structure in (b), indexed by the modified quasi¬ 
momentum —4 k < q < 4 k. The unfolding of the band struc¬ 
ture can be followed from the different coloring of successive 
bands. 


H e ff = ^-fKffW, (3) 

2m 

V e s(x) = cos (2Nkx)a x , U eS = ^ , (4) 

K(t) = —^ sin(2/ca; — f lt)a z H— sm(Nflt)a x . (5) 
nil Nnil 


V(x, t ) = Vl cos(2fcx — fit)a, + Vb cos(Nfit)<j x . (1) 

This potential satisfies V(x + d/N,t + T/N ) = V(x,t), 
with d = 7r/fc, thus, it can be viewed as a continuous 
version of the stroboscopic scheme. Understanding the 
physical effects of the potential (jT]) falls within the de¬ 
scription of time-periodic Hamiltonian systems jl3H16| . 
Following Ref. ca, we describe the dynamics of an atom 
between the times t ; and t[ as 

U(ti -4 t { ) = e - iK< dt) e ~i<dt-ti)H eSe iK(ti), 

where we introduce a time-independent, effective Hamil¬ 
tonian H 0 g and a time-periodic kick operator K{t). The 
three operators in © describe, from right to left, the 
role of the initial phase of the Hamiltonian at time t\, the 
evolution from t; to t[ according to a stationary Hamil¬ 
tonian, and the micro-motion related to the final phase 
of the Hamiltonian at time tf. 

The expressions for the effective Hamiltonian 
H a ff and kick operator K(t) can be calculated 
through a perturbative expansion in powers of 
1/fi, see Refs. 0 R51 . To lowest-order, this yields 


The effective potential Q, which describes a periodic 
potential of depth U e g and spatial period d e g = d/N , was 
derived under the assumption that N is an even integer 
(a similar potential is found for N odd). The expression 
(4) has been obtained based on a Born-Oppenheimer 
approximation, in which the kinetic energy term is 
neglected and one calculates the effective potential for 
a given position x, treating internal degrees of freedom 
(aj operators) quantum-mechanically. One finds that 
the terms neglected here are smaller than those given 
in Q, by a factor E r /(hfi) <C 1. Furthermore, this 
approximation is validated by a direct comparison 
with the full quantum treatment (see below). In the 
Supplementary material we show that the perturbative 
expansion can be resumed, with respect to either the 
variable Vl/( hfi) or Vb /(hfi). There we also discuss the 
generalization of this modulation scheme to an arbitrary 
spin F, through the substitution a u —> 2 F u [2D] , 

In order to test the validity of the effective Hamil¬ 
tonian ©, we performed a numerical study of the full 
time-periodic Hamiltonian using the Floquet formalism. 
Since the Hamiltonian H is invariant under the sym- 
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metries T x '■ x —>■ x + d and Tt '■ t —> t + T, we 
look for eigenstates written as Bloch-Floquet wave func¬ 
tions = e^ qx ~ ut ^u qtU (x,t), where u q ^{x,t) is 

d-periodic in x and T-periodic in t |23][25|. Eigenstates 
are labelled by their quasi-momentum —k<q<k and 
quasi-energy 0 < Hoj < HQ. An example of the band 
structure calculated numerically for N = 4 is plotted in 
Fig#. The band structure exhibits gap openings once 
every four bands, at the momenta Nkp, where p £ Z*, 
as expected for a lattice of spacing d/N. 

The band structure can be unfolded, making use of the 
additional symmetry T* : x —► x + d/N,t —> t + T/N. 
As explained in the Supplementary Material, eigenstates 
associated with the symmetries T x , Tt and T* can be 
written as ipq }U (x, t) = e l( - qx ~ ult ^v qtUJ (x, t ), where v qjU (x, t) 
is d/IV-periodic in x and 27r-periodic in (kx — ttt) [2D]. We 
show the band structure calculated within this formalism 
in Fig. which is very close to that expected for a 
lattice of spacing d/N and depth U e g ~ 10.9 E® ff [26] . 

The practical relevance of the short-spacing lattice de¬ 
scribed above is based on the ability to load atoms into 
the ground state of the effective potential #. The anal¬ 
ysis of this loading protocol requires special care, as the 
effective-Hamiltonian approach inherent to Eq. # as¬ 
sumes a constant lattice depth m ■ In fact, we find that 
the concept of the effective Hamiltonian can be modified 
so as to describe the time-evolution under a ramp of the 
moving-lattice depth Vl, see Ref. [[22] • We simulate the 
lattice loading from a numerical calculation of the full 
dynamics of an atomic wave packet under the action of 
the potential #. Starting from a gaussian wave packet, 
spin-polarized along x, we solve the Schrodinger equa¬ 
tion, discretized in space and time, with a lattice depth 
Vl slowly ramped up for a duration f ramp . As shown 
in Fig. [3^,, a ramp duration t ramp = 20 h/E° s leads to a 
state with strong spatial modulations of spacing d/N, as 
expected for a wavepacket prepared in the lowest band 
of the effective lattice The calculated population 

in the effective lowest band is 93%, close to the value 
expected with standard optical lattices for such a ramp 
duration. In the Supplementary Information we analyze 
the momentum distribution, which corresponds to the 
one expected for the ground state of the effective lattice, 
slightly modified by the micro-motion [2D] . 

The system description as an effective d/N lattice is 
also supported by a numerical simulation of Bloch os¬ 
cillations. We calculate the action of a linear potential 
—Fx applied to the state obtained after the lattice load¬ 
ing. As shown in Fig. 0,c, the wave packet undergoes 
Bloch oscillations, revealed as real-space oscillations of 
its center of mass. Both the amplitude and period of this 
oscillation agree well with those expected for an effective 
lattice of period d/N and depth C7 e flf inferred from band 
structure calculations. 

The potential V (x, t ) written in 0 corresponds to 
the sum of a time-modulated magnetic field and a spin- 
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FIG. 3. (a) Atomic density of a wave-packet loaded into a 
dynamic lattice of spacing d e g = d/4. We start from a gaus¬ 
sian wave-packet, spin-polarized along x, of wave function 
ip(x,t = 0) = exp[—cc 2 /(2cr 2 )], with a ~ 1.4 d (red dashed 
line). The lattice depth Vl is slowly ramped up for a du¬ 
ration tramp = 20 h/E^ s from VL w 0 to Vl = Vl, and lat¬ 
tice parameters N = 4, V® = Vb = Kl = 200F/ ff . The 
atom density after loading is spatially modulated, with a pe¬ 
riod d/4 (blue line), (b) Evolution of the density distribution 
during Bloch oscillations, calculated for the dynamic lattice 
parameters of (a), and for a force F = W e ff/(8d e ff), where 
W e ff ~ 0.06 is the expected bandwidth of the lowest band 
for U e s = 10.9 . (c) Evolution of the center-of-mass posi¬ 

tion during Bloch oscillations, calculated for a standard op¬ 
tical lattice of depth U e g = 10.9 (red dashed line), and 
for the dynamic optical lattice (blue line). The time and spa¬ 
tial coordinates are plotted in units of the ideal Bloch period 
tb = 2Nhk/F and amplitude W e g/F. 


dependent optical lattice moving at the velocity ui a tt = 
H/(2fc). In the above discussion we considered the ef¬ 
fect of this potential as an effective static optical lattice. 
An alternative view is obtained in the frame of reference 
moving at the velocity v = fiatt, where the potential 
V (; x’ = x — vt, t) consists of the sum of a modulated mag¬ 
netic field and a very deep static lattice Vl cos(2kx')a z , 
with Vl ~ hfl U e g. Both points of view can be rec¬ 
onciled by a proper interpretation of the band structure, 
as illustrated for the case N = 2 in Fig. [4] Among the 
eigenenergies u>(q) calculated numerically in the labora¬ 
tory frame v = 0, we identify the Bloch bands corre¬ 
sponding to a static effective lattice of spacing d e g. The 
eigenenergies ui'(q) corresponding to a frame of reference 
moving at a velocity v can be deduced from those in the 
laboratory frame using the relation u/ = lo — qv/h. In the 
frame moving at v = fiatt, we observe Bloch bands cor¬ 
responding to a very deep static optical lattice of period 
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FIG. 4. Band structure corresponding to the dynamic lattice 
for the parameters .V — 2. V'i. = V / b ftfi = 10 E^ s . The 
panels correspond to different frames of reference, of veloc¬ 
ity v = 0 (left) and v = viatt = fl/(2 k) (right). The blue 
points correspond to the band structure of an optical lattice 
of spacing d e ft = d/N and depth U e ft — at rest in the 

laboratory frame. The red dots correspond to the band struc¬ 
ture of an optical lattice of spacing d and depth U ~ 74 E r 
(~ 9 U e fi), at rest in the frame of velocity v = vi a tt |27l . 
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FIG. 5. (a) Momentum-space representation of the effective 
couplings in Eq. [6] illustrated as arrows of length 2 Nk, ori¬ 
ented along the unit vectors ±ei (i = 1, 2, 3), and proportional 
to Pauli matrices. Quantum states are represented in the ba¬ 
sis {|+ 2 ) (filled dots), |— z ) (circles)}, (b) Phase accumulated 
around a triangular subcell of the fc-space lattice. Due to 
the internal-state degree of freedom, the unit cell of the lat¬ 
tice is formed by four triangular subcells. The same phase 
of 4> = 7r/2 is found to be accumulated around all subcells, 
indicating that the lowest energy band is associated with a 
non-trivial Chern number vch = 1 [28]. 


d. 


We now consider a 2D extension of our scheme. The 


time-dependent part of the Hamiltonian is taken as 

V (r, t) = Vl cos(2/cei • r — VR t)a z + Vb cos(NQit)d x 
+ Vl cos(2fce 2 • r — fl 2 t)a- x + Vb cos(NV, 2 t)d y 
+ Vl cos(2fce 3 • r — £l 3 t)d y + Vb cos (NQ 3 t)a z , 

where the unit vectors e 123 have directions as repre¬ 
sented in Fig. [5] For a suitable choice of the frequencies 
^ 1 , 2,3 [29], eac h li ne of the equation above can be treated 
individually, which results in an effective potential of the 
form [30j 


T4ff(r) ~ — V-[cos(2Nke 1 ■ r)a x + cos(2 Nke 2 ■ r)a y 

+ cos(2Nke 3 ■ r )a z ], (6) 

where N is taken to be an even integer. These couplings 
are illustrated in quasi-monrentum space in Fig.[5ji. Fol¬ 
lowing Ref. PS], the topological Chern number associ¬ 
ated with the lowest energy band can be readily obtained 
from these couplings. Indeed, the Chern number mea¬ 
sures the flux of the Berry curvature f2(q) over the entire 
(momentum-space) unit cell: 

^Ch = 7 T- f fl(q)d 2 q, (7) 

^ J unit cell 

which can be directly evaluated by calculating the phases 
accumulated by a state as it performs a loop around 
the triangular subcells [25]. For the effective lattice de¬ 
scribed by eq. § , each unit cell is constituted of four 
triangular subcells, and we find an accumulated phase of 
7 t/ 2 within each of them (see Fig.jsjr). In this configu¬ 
ration, the Chern number of the lowest band is given by 
^Ch = (1/27 t) x 4 x (tt/ 2) = 1. Generally the reasoning 
above is valid only in the weak-binding regime; however, 
for the coupling ([6]), ^ch is unchanged for all values of 
U e g. Note that the size of the unit cell in real space 
scales as 1 /N 2 \ we thus expect the flux density to be in¬ 
creased by a factor of N 2 compared to standard optical 
lattices. 

The method discussed above is based on applying 
strong spin-dependent optical lattices, for which Lan¬ 
thanide atoms would be most suited for a practical im¬ 
plementation. Indeed, the optical lattices could be gener¬ 
ated using laser light close to a narrow optical transition, 
which would lead to deep, spin-dependent lattices with 
negligible Rayleigh scattering effects (for Dy atoms, one 
can achieve ratios ftFgcattering/Vi, ~ 10 -7 ) [311433] . 

In conclusion, we introduced a novel scheme to en¬ 
gineer spatially periodic atom traps of sub-wavelength 
spacing, based on the application of spin-dependent op¬ 
tical lattices. A natural extension of this work would 
be to include interactions between atoms in the effective 
lattice description, and to understand whether micro¬ 
motion plays a significant role in scattering properties 
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[MHMj . This aspect will play a central role for investi¬ 
gating quantum many-body physics with short-spacing 
lattices. 
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S.I. EFFECTIVE POTENTIAL CREATED BY 
THE STROBOSCOPIC SCHEME 

In this Section, we present the calculation of the ef¬ 
fective lattice potential produced via the stroboscopic 
scheme illustrated in Fig. 1 (main text). We start from 
a periodic potential V ( x ) of period d, which can be de¬ 
composed in Fourier series as 


and discuss further the relevance of this approximation 
for typical lattice parameters. 

The Born-Oppenheimer approximation consists in ne¬ 
glecting the kinetic energy for the calculation of the effec¬ 
tive Hamiltonian. The position x is considered as a fixed 
parameter, while internal degrees of freedom are treated 
quantum mechanically. Using the formalism of ref. [SI], 
we decompose the Hamiltonian in Fourier series as 


V(x) = J2 V r 

pel 


gi27vpx / d 


The method consists in shifting the potential V(x) of 
the distance d/N, after each time interval T/N. For N 
integer, this leads to a time-periodic potential V(x, t) 
of time period T. For a sufficiently short period T, 
the atomic motion is governed by the effective potential 
V!eff(x), equal to the time average of V(x,t): 


U e ff(a:) 


T 


V(x,t)dt 


N 


-^(x + jd/ N) 


3 =1 


N 


pgZ j =1 

V p e l2 * px/d . 

p multiple of N 


i2irpj /N 


It is then apparent that the effective potential V e s(x) is 
periodic, of period d e ff = d/N. 


H(t) =H 0 + J2 V U) e ijnt + 

3=1 

V {1) = U ( ~ 1)t = V h e~ 2ikx a z / 2, 
V (jv) = U (-JV)t = V B a x /2, 
yb) = 0 otherwise. 


The effective Hamiltonian can be expanded as a series 
in 1/17. At lowest order, it reads 


H eS = 
V eS (x) = 


2 m 


Wff (*^) 7 




N\{Ml) N L 


V U) V (-N) 


h.c. 


with VW occurring N times. One calculates 
V N V B e~ i2Nkx 

Ves{x) = 2N\(2hUl) N [CT - ' •' ’ [<Tz ’ + h ' C ' 

_ U eff J cos(2Nkx)a x , N even, 

2 \ sm(2Nkx)a y , N odd, 


S.II. CALCULATION OF THE EFFECTIVE 
POTENTIAL 

In this section we give details about the establishment 
of the effective Hamiltonian (eqs. (3), (4) in the main 
text). We first provide a simple description of the effec¬ 
tive Hamiltonian in a Born-Oppenheimer approximation, 


* sylvain.nascimbene@lkb.ens.fr 


where U e s = 2V^{V\,/hUL) N /N\. 

In the main text we consider the case of a (pseudo)- 
spin-1/2 for simplicity. The scheme can directly be ex¬ 
tended to an arbitrary spin F. Assuming a coupling 
V(x,t) = 2Vl cos(2fcx — Q,t)F z + 2Vb cos(NClt)F x , we 
make use of the general commutation algebra of F u 
operators and obtain an effective potential V e g[x) = 
U e g cos(2Nkx)F x (N even). 

We now consider the validity of the Born-Oppenheimer 
approximation, which consists in neglecting the non com¬ 
mutation of position and momentum. The latter plays a 
role at the order N + 1 of the perturbative expansion 
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v L [fin] 



v B [fin] 


FIG. SI. Depth U e a of the effective lattice (in absolute value), 
calculated for arbitrary values of Vl/(KI) (a) or Ve/(fin) (b) 
in the case IV = 4. The dashed lines correspond to the lowest- 
order perturbation result (4), and the solid lines to the resum¬ 
mation results (S.l) and (S.2). 


in 1/n. We performed the perturbative expansion up to 
the order N + 1 for the case N = 1, and obtained the 
expression 

H e g =^— + ^ sin [2kx)a y 
2 m 2 

+ (j|) E ’ + ^f^n <!> si ”< 2fa ') + Sin(2tx)p) 

While the first term of order N +1 only represents an en¬ 
ergy offset, the second term can be viewed as a spin-orbit 
coupling which may affect the atom dynamics. However, 
for typical momenta p on the order of the lattice momen¬ 
tum k , the amplitude of this term is smaller than the 
effective lattice depth U e s by a factor ~ E T /(Ml), which 
is typically very small for the examples considered in this 
Letter. 


S.III. RESUMMATION OF 
THE PERTURBATIVE EXPANSION OF H eS 


The effective potential U e ff can be calculated as a se¬ 
ries expansion in powers of the (potentially small) dimen¬ 
sionless parameters Vl/(MY) and V B /(hfl), in the high- 
frequency limit Q —► oo. In the main text, we provide its 


expression in Eq. (4), which corresponds to the lowest- 
order term. We note that this derivation, which is based 
on the general formula of Ref. [SI], was obtained by ne¬ 
glecting the non-commutativity of the position and mo¬ 
mentum operators; indeed, we verified that the momen¬ 
tum operator is irrelevant in the derivation of the effective 
potential, which essentially relies on the spin-dependent 
time-modulated components of the Hamiltonian. Thus, 
in the following of this Section, which aims to derive the 
effective potential in the strong-modulation regime, we 
explicitly neglect any effects associated with the kinetic 
energy term of the full Hamiltonian. 

In this Section, we first derive the expression for the 
effective potential V e s, in the case where Vl/(M l) is al¬ 
lowed to take arbitrary large values (still assuming that 
V B Vl , fiO). Following Refs. [S2, S3], we perform a 
unitary transformation 


1 ^') = R(t) W ), 


R(t) = exp 




which removes the diverging term ~ Vl ~ hfl from the 
time-dependent potential V(x,t) in Eq. 1 (main text). 
This leads to a novel time-dependent potential 

V'(x,t) = R{t)V(x,t)(t)R\t) + ihd t R(t)R}{t) 

= R(t) [V B cos (NSU)<t x ] R*(t). 

Making use of the identity a x e lirT * = cos(2j)a x + 

sin( 27 )<T y , we obtain the expression 


V'(x,t) = V B cos(Nflt) 


cos ( ^y sin (kx — f It) ) a x 
2y l 

+ sin ( sin(fca: — fit) ) a y 


In the large-frequency limit fl —> oo, the atom dynam¬ 
ics can be described by an effective stationary potential, 
given by [S2, S3] 


V e g{x) 


Jo v '^ x,t ^ dt 



Vb cos(2Nkx)a x , 


(S.l) 


assuming N even, and where Jjv is a Bessel function of 
the first kind. This effective potential corresponds to a 
spin-dependent optical lattice of spacing d/N and depth 
l-\,< = 2 ./.V ® Vb. 

A similar resummation with respect to V B / (Ml) can 
also be derived. Here, we make use of the Floquet repre¬ 
sentation of time-periodic Hamiltonians. We first write 
the exact eigenstates of the coupling ^ cos (Nflt)cr x , 
which read 


I n, s : 


■)- 52 Jp 


pe z 


(2s x V b \ 
V Ml ) 


I n+pN, s x ), 
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where n denotes the Floquet quantum number, and s x 
is the spin projection along x. The energy of the state 
||n, s x ) is equal to nMl. The effect of the coupling 
Vl cos(kx — flt)a x can be understood using perturbation 
theory in the degenerate subspace ||n, ±), which must be 
performed at order N. We obtain the expression 


VrffO) 


Uefi 

2 


cos(2 Nkx)a x , 


u eS = 4 nn 



E «^ = - 1 n^ 1 E;=i( 1 + Ar ft) 

(S.2) 


We plot in Fig. SI the lattice depth U e & given by the re¬ 
summation formulas in Eqs. (S.1)-(S.2) discussed above. 
We checked that the formulas (S.l) and (S.2) account 
well for the numerical results obtained via direct di- 
agonalization of the Bloch-Floquet equations (see Sec¬ 
tions.IV). 


S.IV. EXPRESSION FOR 
THE BLOCH-FLOQUET HAMILTONIAN 

The modulated potential (1) is invariant under the 
space and time translational symmetries 7~ x , 7* and T*, 
which all commute with each other. The eigenstates of 
the Hamiltonian can thus be written as eigenstates of 
those symmetries, which can be expressed as 

tpq, u (x,t) c jtl e ilikx ~ nt) e ijNkx , 

3,iez 


where — Nk < q < Nk and 0 < ui < Q. The spinor 
coefficients Cjj are determined by the equations 


_h 2 [q + {l + Nj)k} 2 ^ 

Ti\Oj l&ljCjj — ^ CjjZ 




+ -^ cr x{ c j,l+ 1 + c j,l-1 ) 

V B 

+ —^(Cj+i^-JV + Cj-l,j+jv)- 


(a) 3 

'J 2 

7 ? 1 


(b) 2 
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FIG. S2. (a) Momentum distribution associated with the 

ground state of the effective lattice with spacing d/4 and 
depth U e s = 10.9 E^ s . (b) Momentum distribution of the 
state in (a), taking into account the micro-motion expected 
for the dynamic lattice parameters, according to Eq. S.3. The 
micro-motion leads to a more complex structure compared 
to (a), periodically evolving in time. The lattice parameters 
correspond to the ones of Fig. 2 in the main text. 


t = 0 

1 1 1 | i i i | 


1 1 1 1 1 1 i i 

t = 0.25 T/4 

1 1 1 | 1 1 1 | 


t = 0.5 T/4 

1 1 1 | 1 1 1 | 


t = 0.75 T/4 

1 1 1 | i i i | 

Ly,V— r— ^ , T . T + . 


-8 -4 0 4 8 

q [2k] 



the family of states of momentum multiple of 2Nk (see 
Fig. S2a). 

The actual state created using time-modulated lattices 
is expected to be modified by the micro-motion, as 

I m) = e~ lK{t) IV’o) , (S.3) 


The numerical data represented in Fig. 2 (main text) is 
calculated using the above equations, in a truncated basis 

-10 <j,l< 10. 

S.V. MICRO-MOTION EFFECTS 
IN THE MOMENTUM DISTRIBUTION 


where the expression for the kick operator K(t) is given 
in the main text [Eq. (5)]. The latter leads to additional 
diffraction peaks at all momenta multiple of 2k, whose 
amplitude vary periodically in time, with a period T/N 
(see Fig. S2b). This shows that Bragg diffraction does 
not give a direct information on the ground state of the 
effective lattice. 


In this section, we analyze how the micro-motion as¬ 
sociated with the time-modulation in Eq. (1) affects the 
momentum distribution of atoms prepared in the effec¬ 
tive potential V e g of spatial period d/N [Eq. (4)]. Specif¬ 
ically, we consider an atom prepared in the ground state 
of the effective potential. This state can be expanded on 


S.VI. EFFECTIVE HAMILTONIAN 
DURING LATTICE LOADING 

In this Section, we analyze the adiabatic preparation of 
the ground state associated with the effective potential 
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FIG. S3, (a) Evolution of the amplitude of the moving optical lattice at position x = 0, during and after the lattice ramp of 
duration t ramp . (b) Scheme of the ramp discretization: the time interval 0 < t < tramp is decomposed into A f steps of duration 
At. Within each step the depth Vl is constant, leading to a time-periodic potential. 


V e s(x) of depth U° e . We consider a slow ramp of the 
moving-lattice depth Vl (t) = 0 —> V® during the time 
interval 0 < t < t ramp , such that the effective potential’s 
depth U® s corresponds to the final value Vl (tramp) = V°. 
As the definition (2) of the effective Hamiltonian and 
kick operators assumes a constant lattice depth [SI], we 
expect these notions to be modified during the ramp. It 
is the aim of this Section to show how the adiabatic ramp 
can still be captured by an effective-Hamiltonian picture. 

To analyze this situation, we decompose the ramp into 
A f steps, and we assume that the time interval At = 
tramp/A/ is short enough, such that Vl can be considered 
to remain constant within each step. More precisely, we 
assume that the lattice depth is equal to VL(jAt) during 
the step jAt<t< (j + 1)A t. We then apply the effective- 
Hamiltonian formalism of Ref. [SI] within each time-step, 
and write the full time-evolution operator as 

o 

fAamp — 11 Uj , 

j—M-1 

JJ. — e -i K 0 [VL (j At)] e -iH e ti[V h {jM)]M/n e iK 0 [VL O'At)] 

In the latter expression, and for the sake of simplicity, 
we assumed that At was a multiple of the modulation 
period, so that the kick operators at the beginning and 
at the end of each step only depend on the value of Vl ( 
in fact, they correspond to the kick operator at the time 


t = 0, hence the notation Kq). 

Assuming At short enough, we write 

e iK 0 [VMj+l)At)] e -iK 0 [V^U&t)] _ e iAi(dV L /dt)dJWdV L 

leading to 

Ur amp = e _iif(tramp) T jexp J (t)dt / hj | , 


where T denotes time-ordering, and where one intro¬ 
duced the slowly varying Hamiltonian 


^ cos(2Nkx)a x 


VlW 

sm(2kx)a z . 

(S.4) 


We now estimate a criterion for identifying the adia¬ 
batic regime of the lattice loading. Describing the lattice 
loading solely with the first term of eq. (S.4) would lead to 
the standard adiabaticity criterion f ramp HU 9 S /(E^ S ) 2 
[S4]. We expect the second term in (S.4) to drive non- 
adiabatic transitions for Vl > ClE r . Adiabatic lattice 
loading thus requires the additional constraint f ram p ^ 
(Vl /E t )/£1. As we choose Vl < hfl, this constraint 
should not be the most restrictive. 


[SI] N. Goldman and J. Dalibard, Phys. Rev. X 4, 031027 [S2] P. Hauke, O. Tieleman, A. Celi, C. Olschlager, J. Si- 

(2014). 
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